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A microscopically based Hamiltonian of the generalized 
t-t'-J model is presented. Two types of the additional t' - 
terms are discussed. The numerical range of the ampli- 
tudes corresponding to the additional i'-terms for the real 
CuC>2 planes is derived from the three-band model calcula- 
tions. Using the variational spin-polaron approach the single- 
carrier dispersions in the generalized t-t'-J model are calcu- 
lated both for the hole- and electron-doped systems. The 
hole and electron band minima are found to be at points 
(±7r/2, ±7r/2), (0, ±tt), and (±tt,0), respectively. The band 



minima shifts |A 



(0,7r)-(7r/2,ir/2) 



are not small (~ J). The 
bandwidths for both cases of doping are found to be 1.5 — 4.0 
times larger than those in the t-J model. 
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I. INTRODUCTION 

There is a general agreement that the three-band Hub- 
bard model is an appropriate basis for consideration of 
the spin and charge excitations in the Cu02 planes of 
high-temperature superconductors One of the in- 

teresting problems for these systems is the hole energy 
spectrum. One-hole energy calculation in the framework 
of the above mentioned model has been done by Bara- 
banov, Kuzian, Maksimov, and Uimin who used 

the variational approach. This method seems to be too 
complicated since the characteristic energy scale of the 
three-band model is a few electron-volts while that of the 
hole energy spectrum is several tenths of an electron-volt. 
Therefore, it is natural to obtain first the low-energy limit 
of the three-band model |]6|-|l2| and then to investigate 
the spin and charge degrees of freedom in the framework 
of this effective model [|l3| {[7) . 

It is widely accepted by now that the simple model 
which contains in itself the interacting spin and charge 
degrees of freedom is the so called t-J model 



H t -j = t 4, a cj, a + SiSj, 

(ij) 



(1) 



where (ij) denotes the nearest neighbor sites, Si is the 
local spin operator, c\ a {ci t0 ) is the constrained fermion 
creation (annihilation) operator. The explicit form of 
this constraint depends on the type of doping. 



The considerable attention that this model has re- 
ceived originates partly from the fact that it can be de- 
rived from the above mentioned three-band or more gen- 
eral d — p models which take into account the detailed 
electronic structure of the copper oxides. Recently it 
has been shown pi that this derivation can be done 
quantitatively for the real Cu02 planes by the use of an 
additional calculation of some experimentally observable 
values. 

The behavior of a single quasiparticle in the t-J model 
has been studied intensively using both analytical and 
numerical techniques |]l4] , |l5[^9pc| l . These investigations 
have clearly shown that the naive tight-binding picture 
is completely inadequate for the carrier motion on the 
antiferromagnetic background. Namely, for the realistic 
t/J ~ 2 — 3 the bandwidth is of the order of 2J, not 
Wq — 8t. This strong (of the order of 10 times) band- 
width suppression results from the distortion of the spin 
background by the carrier hopping from one sublattice 
to the other. In this situation, when the characteris- 
tic energy scale becomes J, even small (compared to t) 
single-sublattice hopping t' can be the key parameter for 
the subtle details of the energy spectrum and other fea- 
tures of the doped systems. This is clear since the single- 
sublattice motion does not disturb the spin background. 
Hence, a careful analysis of the low-energy single band 
limit of the realistic d — p models should be done keeping 
all essential terms over the t'/J (not t' /t) parameter. 

The simplest form of the additional i'-terms can be 
written as 



Ht , = t> £ 



(2) 



emphasizing that the nonzero 0-0 hopping provides 
large enough transfer amplitude to the next-nearest 
((27)2) neighbor |2(],|2l]]. This term (||) alone has been 
found to be responsible for the Cu02 plane electron-hole 
asymmetry and to be useful for the interpretation of the 
recent angle-resolved photoemission experiments [ p2~| . 

The goal of the present paper is to exhibit the gen- 
eral form of the additional i'-terms which one can get 
from the three-band Hubbard or the other first-principle 
models, to determine numerical ranges of all essential 
parameters for the real CUO2 planes, and to calculate 
a single-particle dispersion using the spin-polaron ideas. 
Namely, we will demonstrate that two types of the ad- 
ditional i'-terms naturally originate in the single-band 
model: (i) the terms arising from the 0-0 hopping, and 
(ii) the second-order high-energy channels terms [E3--E5) . 
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The importance of both of them will be clearly shown. 
The first ones are the key parameters for the band min- 
ima shifts, while the second ones lead to the growing of 
the bandwidth. A quite similar model was investigated 
in the mean-field approximation for a special range of 
parameters in the work by Onufrieva et al [^6| . This ap- 
proximation neglects some essential features of the hole 
(electron)-spin interaction and thus has only a qualitative 
character. 

Single-hole energy calculation in the limiting case of 
the three-band model using the self-consistent Born ap- 
proximation has been done in the work by Starykh et al 
p7[ . In spite of neglecting some features of the problem, 
the results of this approach are in a qualitative agreement 
with the present paper. 

The quantitative reduction of the three-band model 
developed in Ref. Jji| for the real CuC>2 system provides 
the method of the derivation of a set of the realistic ranges 
for all additional i'-terms. We will show that the band 
minima shifts as well as the bandwidths are not small for 
these ranges of parameters. Also, the density of states 
(DOS) characteristic features will be discussed. 

The paper is organized as follows. In Sec. II we 
give the general form of the corrections to the Hamilto- 
nian ([!]), discuss them, and show the possible parameters 
ranges. In Sec. Ill we discuss our results for the spin- 
polaron dispersions and for the DOS features. Finally, in 
Sec. IV, we draw conclusions. 



II. REALISTIC LOW-ENERGY MODEL 

The most accurate form of the effective f'-terms that 
follows from the three-band Hubbard model is p8| , p8[ : 

E ia^ + t" £ <<A* 

+ t N E A^o^+t 3 ^pa^pSi, (3) 



(ilj),a 



(ilj),a/3 



where (27)2(3) denotes the second (third) next-nearest 
neighbor sites, (ilj) denotes the three nearest neighbor 
sites ({il), (jl)), Ni is the number of fermions operator, 
S/ = a <J aj3 c l,0 is the local spin operator, er is the Pauli 
matrix, a = —a. The differences in the explicit form of 
the ct (c,- a ) operators for the hole- and electron-doped 
systems as well as the signs of t',t",t N ,t s will be dis- 
cussed later. Now we will consider the origin of the t', t" 
and t , t s terms. 

Let us begin from the last ones. The second-order per- 
turbation treatment of the usual Hubbard model near 
half filling provides the effective superexchange interac- 
tion (Eq. ph) and the so called 'three-site spin-dependent 
hopping' p4,py pl| • Their importance for both the spec- 
trum of the charge excitation and the hole-hole interac- 
tion was noted p9| . These terms in Eq. (||), manifestly, 
have the rotationaly invariant form. In the case of the 



usual Hubbard model mapping to the t-t'-J model, the 
expressions for the t , t terms have the simplest form: 



t N = = 



1 



2U 



J 
8' 



(4) 



It should be noted that this expression (|J) is valid for 
both types of doping, and that the signs of t N and t s are 
the same for the electron and hole due to the second-order 
nature of these terms. What would one expect from the 
three- or more-band Hubbard model mapping? 

Much more high-energy channels (triplet, etc.) are 
opened for the virtual second-order transitions in com- 
parison with the usual Hubbard singlet only case 



This leads to 

N,S _j_ +N,S 



two effects: 



(i) t N ± -\t b , and (Ti) 
where e and h denote the electron and hole 



hopping integrals, respectively. The exact expressions for 



(t 



e,h 



and (t ) e ,h were derived in Refs. Ql 



In spite 



of the contribution of the highest states, the lowest sin- 
glet is the most important. The main features of these 
terms in the effective Hamiltonian (^) remain unchanged, 
i.e.: (i) sign(t N ) = —sign(t s ) = +1 both for the electron 
and hole, and (ii) \t N \ ~ \t s \ ~ J/4, so that they are not 
negligible. 

In addition, one can point out that at the mean-field 
level, N and S operators should be replaced with their 
averages (N) and (S), which effectively leads to: 



H 



N,S 
MF 



? E 

<*J>2,< 



St 7,. 



C i,a C j,a + t 



E 

<»J>3>< 



(5) 



with t' = 2t" = 2((N)t N + (S)t s ), where the additional 
factor 2 for the diagonal hopping integral t' as compared 
to the 'oversite' t" arises from the two possible ways for 
the virtual processes on the square lattice. Thus, the 
second-order terms in the Hamiltonian (||) can be approx- 
imately considered as the renormalizations of the 'bare' 
t',t" amplitudes (||). As was noted, these renormaliza- 
tions do not change the signs under changing of the dop- 
ing type. 

Now, return to the first two terms of the Hamiltonian 
(^). They are the first-order terms arising in the first- 
principle models for the Cu02 plane from the nonzero 
0-0 hopping. While including of the i'-term is evident 
Jl5| , p0 21 1, including of the i"-term requires an additional 
explanation. Physically, only the i'-term can arise for 
the next-nearest neighbor CUO4-CUO4 local states due 
to the tjm matrix element (t]]3^]. However, as was shown 
HUj^l^l the correct state of the oxygen low-energy de- 
grees of freedom are the ortogonalized Wannier states. 
Consequently, some 'unphysicaP transition amplitudes to 
the more distant neighbors arise. These amplitudes fall 
rapidly with distance, and only the t', i"-terms should be 
kept in the low-energy model Q . 

As was recently proposed in many works, the absolute 
sign of the £'-terms for the electron- and hole-doped Cu02 
systems could be the source of the strong electron-hole 
asymmetry in the magnetic phase diagram 0,p4| . 
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The electron-hole asy mmetry arises naturally in the 
three-band model @,@|U||]. Not only the signs of the 
hopping integrals, but also the orders of the elementary 
processes, that lead to the effective hopping, differ for 
the electron and hole. Consideration of a simple limiting 
case as well as discussion of the electron-hole asymmetry 
and the dependence of the effective hopping integrals on 
the initial parameters are presented in Appendix A. Nu- 
merical results for the ranges of the effective parameters, 
based on the exact formulae of Ref. j2f| and the approach 
developed in Ref. |l8[ , are shown in Table 1. 

Further, we will consider these (Table 1) t',t",t N ,t s 
values and th = 2.5 J, t e = — 3.0J [^2|,@ as the realistic 
ones. 



III. SINGLE-CARRIER ENERGY SPECTRUM 

Before doing the energy spectrum calculation let us 
discuss the sense of the constrained Fermi operators ct(c) 
introduced in Eq. ([[]) . Both the hole and electron single- 
band low-energy Hamiltonians derived from the three- 
band Hubbard model are naturally expressed in terms of 
the Hubbard operators at the site i 

Xf = \ai)(bi\, (6) 

here the states a, b are either spin or singlet (vacancy) 
local states. Since the undoped state of the CuC>2 plane 
is set up by the localized holes in the Cu d 10 and O p G or- 
bitals, the additional low-energy hole forms a singlet with 
the local one, whereas the electron in this background is a 
vacancy. Therefore, it is convenient to rewrite t, t'-parts 
of Eqs. (|), (|) as: 

H hole = y $ j xrxf+ y tfxrXj'fit 

(ij)n,a {ilj),af3 

+ Y, fiXrXf'iVflaSl), (7) 

(Uj),a0 

H el = Y te ij x i ax f + Y t?X? a XfN t 

(ij)n,a (ilj),a/3 

+ Y W*fV*,S,), 

(ilj),a0 

where |a) = | f), | J.) is a local spin state, \s) and |0) are 
singlet and vacancy states, respectively, (ij) n denotes all 
essential neighbor sites. Note, that the signs of the next- 
nearest hopping parameters shown in Table I are related 
to these particular (Eq. (0)) order of the Hubbard op- 
erators. The signs of the nearest neighbor hopping are: 
\v) > ^' < ®' ^ e ^ a ^ er ^ s n0 ^ essen tial, since for 
the antiferromagnetic background the shift of the quasi- 
momentum space by the vector of the reciprocal lattice 
Q — (-7T, 7r) changes the sign(tuj\) but physically changes 
nothing jlq ]. 

We wish to stress that when the kinetic-energy part 
is of the single-band type (as in Eq.(^)), mapping of the 



Hubbard operators onto the constrained Fermi operator 
basis is twofold. 

If one prefers to retain the singlet ideology for the hole- 
doped system, the first term of the Hamiltonian (fjj) can 
be rewritten identically as H |: 

H hole = y tiA-A^rf,-^ ( 8 ) 

and the background is created by the hi Q -hole at ev- 
ery site. n\ _ a = h\ _ a hi_ a is the projection operator 
which project out the vacancy states ('double electron 
occupancy'). This representation for the electron over 
the hole background is: 

H * = - E W-nt-aKvhjAl-^-a), (9) 

(*J>n,a 

where the operators (1 — n,;. Q ) project out the double 
hole occupancy. An alternative way is to consider the 
physical hole as a 'hole' (vacancy) in the upper Hubbard 
band |3^,^], and the physical electron as the 'particle' 
(singlet) in the lower one. This freedom in choosing is 
not connected with the initial orbital structure of the 
CuC>2 plane, but follows from the algebra of the Hubbard 
operators in Eq.(Q). Hence, Eq.(^) for an extra electron 
on the electron background can be written as: 

Hel = Y t V n i-<x e la e 3,'x n j-a- ( 10 ) 

At the first glance, it would seem that the choice of the 
'particle' (Eq.([fO])) or 'vacancy' (Eq.@) language leads 
to the change of the sign of the quasiparticle energy. This 
is not true, since the 'vacancy' energy has the sign of Uj 
reversed compared to the 'particle' energy |T^,^I|. There- 
fore, Eqs. (J9j) and ( |Io| ) lead to the same energy. 

Thus, the explicit form of the constrained operators in 
Eqs. (0), @ is: c itC , = h^ a n\ a {ei^n\^) for the physical 
hole (electron) system. Hereafter, we will work with the 
Hamiltonian Eqs. (Q), (^), constraint from Eqs. (p|), (pl|), 
and parameters from Table I. 

Properties of the single-particle in an antiferromag- 
netic background were studied in detail by many authors 
|,|g|Jl||rj,^-|9| using different approaches. Their re- 
sults coincide at the point that the carriers are strongly 
dressed by the spin waves, i.e. quasiparticles are the 
magnetic polarons of a small radius with a strongly 
anisotropic dispersion law and small enough quasiparticle 
residue. We base our calculations on the results of Ref. 
|l4| . The suggested trial function of an extra particle has 
the simple form p[ : 

+ si Y vk,s4i+6,o ex p(* kr «), 

<5 
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valid also for , after changing f-o-J. and S + <^4> S . 
n Gf (4) means the site n of the spin-up (-down) sub- 
lattice, 5 is the unit vector to the nearest neighbor site. 
The explicit expressions for and /x kj <5 are given in Ap- 
pendix B. From the string picture point of view this sim- 
plest ansatz ([ll]) consists of the 'bare' particle and four 
shortest 'strings' of the lattice constant length. It was 
shown |Il| , f40| -[42|] that the one-hole dispersion as well as 
two-hole contact and long-range interactions, and even 
many-hole properties of the t-J model are quantitatively 
well described using the above ansatz. 

Consideration of the generalized i'-terms (Eq. (j3j) ) does 
not require any changes in the trial function (11). 
Roughly, this is due to the following reasons: (i) the t'- 
terms lead to the motion of the particle over one sublat- 
tice without distortion of the spin background, that en- 
ables 'bare' particle to propagate freely; (ii) the motion 
of the 'dressed' particle (with a string) mainly leads to 
longer strings, whose contribution to the energy is of the 

order of and parametrically small at t ~ 3 J. Hence, 
the expansion of the ansatz makes a small decrease in the 
energy, since the part of the ansatz ( |lT| ) , which would be 
mainly affected by the t'-Hamiltonian (||), is the 'bare' 
one. The last statement will be demonstrated below. 

As was noted, the i'-induced transitions occur in one 
sublattice, therefore a good preliminary consideration of 
the role of each term in Hamiltonian (j|) can be done 
for the Ising background. It allows us to find easily the 
leading contribution to the band minima shifts, band- 
widths, and effective masses. The N and S operators in 
(g) should be replaced with (N) = 1, (S) = (S z ) = ±1/2. 
The eigenenergy of the magnetic polaron (ujh is 

McKfft-j + gtOh/fc) 2 a 

(WW 



difference A 



(0,ir)-(ir/2,7r/2) 



02 is proportional only to 



X = Av 2 (t' + 2t'' 



At 



N 



2t% 



02 =4v 2 (2t" -t>), 

7k = l/2(cos(fc x ) + cos(k y )), 

7 k = l/2(cos(fc x ) - cos(fcy)), 



(12) 



where we omitted all high-order (t' e ff/t) terms. Eo is the 
depth of the band, 0\ , 0i are the inverse masses in the di- 
rections ortogonal and parallel to the magnetic Brillouin 
zone boundary, respectively. Since for the Ising case a 
particle of the pure t-J model is dispersionless |i3|], v 
and \i are constants and at the realistic t ~ (2 — 3) J 
i/| 2 ~i, |/i| 2 ~i. Expression ( |l2] ) for the energy co- 
incides almost exactly with that obtained for the free 
spinless fermions |^6|. The main difference is the pres- 



ence of the weight of the 'bare' particle \i>\ 2 in Eq. (12) 



In the realistic region of t, t' , t" etc Eq is very close to its 
t-J z model value E — —It. 

Two notes should be done. Firstly, the /3i-term jl2] ) 
represents the dispersion which is degenerate along the 
(tt, 0) — (0, 7r) line, whereas the /3 2 -term lifts out this de- 
generacy, placing the minima at the (±7r/2, ±7r/2) (if 
02 > 0) or (±7r,0),(0,±7r) (if 02 < 0) points. The 



the 'first-order' terms t',t" (12) and does not depend 



on the second-order t , t terms. Secondly, one can 
sec from Eq. (jl|) that in the case t N > 0, t s < 
and sign(t') = —sign(t") 0\ and bandwidth W(= 0\ 
or 0i —02) are mostly determined by a few (2t N — t s ). 

The above qualitative calculation for the simple Ising 
background results in conclusions, which remain valid 
for the Neel case. The energy difference between the 
(7r/2,7r/2) and (w, 0) points, which is the crucial value 
for the various calculations, has different signs for the 
hole- and electron-doped systems and is not small for 
both cases. It means that if the single-particle picture is 
valid for the finite doping regime, the quasiparticle Fermi- 
surface will be located (up to the high enough doping 
level) near the points (±7r/2, ±7r/2) and (±7r,0), (0, ±7t) 
for the hole- and electron-doped systems, respectively. 

Also, it can be shown from Eq. (12) that, at least for 
the upper limit of the t' ,t" ,t s ,t N hopping parameters 
(Table I), the bandwidths are large enough. Thus, one 
would expect decrease of the DOS compared to the pure 
t-J model one. 

Finally, the role of the 'bare' particle for the consid- 
ered simple case is crucial, since the 'dressed' part can- 
not propagate freely in the absence of the spin fluctu- 
ations. Therefore, the weight of the bare part in the 
trial wave function is the subject of the prime interest. 
As was noted n bare = v 2 ~ 0.5, that is larger than in 



other works |3jj-|3g] where n bare = v 2 ~ 0.35 - 0.41. 
For n lstrins = 4|^| 2 ~ 0.42 the agreement is better: 
filstnng ^ q.44. This discrepancy is due to the transfer of 
the weight from the rest of the exact particle wave func- 
tion (with the infinite number of strings) to the weight of 
the 'bare' particle in the approximate ansatz (11). Thus, 
our main approximation lies not in the 'shortness' of the 
ansatz, but in the overestimation of the bare particle 
weight. This problem can be overcome by the simple 
renormalization of v 2 to v 2 . 

Let us discuss now what changes of the above results 
can be expected for the Neel background. The main 
changes arise from the fact that due to the spin fluctua- 
tions the pure t-J model particle is given the possibility 
to propagate. Evidently, this will provide additions to 
0i,02 and the bandwidth W. Also, the 'dressed' part will 
lead to some coherent transitions. Due to the more com- 
plex structure of the matrix elements of the Hamiltonian 
(0)) @ f° r tne Nfel background, the simple tight-binding 
relations between 0x,0i and W no longer hold. 

The pure t-J dispersion law [^4pl{| is: 



e*" J (k) = 2J - yj (0.66 J) 2 + 4.56i 2 — 2.8i 2 7 k 



+0.01|t|( 7k ) 5 



(13) 



E, 



t-j 



t-j, 



where 0[~ J ~ +0.65|t| and 0l~ J ~ +0.01|t|. Due to 
the k-dependence of the trial wave function components 
^kj/^k and the nonzero matrix elements for the string 
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components, the first order contribution to the energy 
from the t'-Hamiltonian (^|) contains extra terms with 
the highest powers of 7^: 



e(k) = e *- J (k) 
<5e 4 '(k) = SE Q + S0 1 ^ 

+ ^k7k + S k7 ^(7 k ) 2 



(14) 



^2,k(7 k ) 2 



where the 5/3i,k 5 <5/?2,k, ^k and £?k coefficients weakly de- 
pend on k through the |i^k| 2 and | /u-k | 2 quantities (see 
Appendix B). As was discussed earlier, the highest-order 
corrections to the energy from Eq. (^) are of the order of 

, and for the realistic t/ J ~ (2 — 3) and t! from Table 
I can be neglected (~ j^). Since |^k| 2 is almost constant, 
the renormalization to the 'normal' \u\ 2 can be accom- 
plished by the simple replacing 1^^/2,^/2 1 2 — l^.ol 2 — 
Ko| 2 ^ |i>| 2 ~ 0.35. 

Figures present our results for the hole and electron 
dispersions. Contour plot in Fig. [l] shows the character- 
istic feature of the hole dispersion in the whole Brillouin 
zone. The magnetic Brillouin zone boundary as well as 
the r(0,0) -> M(7r/2,7r/2) -> Z(tt,tt) -> V(0,tt) -» T 
and F-tM-» X(ir,0) directions are shown. Figs. ^||| 
show the dispersions along the r — ► M — > Z — > Y —> T, 
and Y — > M — > A lines. Fig. || demonstrates the charac- 
teristic t-t'-J hole dispersion (solid line) for the average 
parameters from Table I, pure £-J model hole disper- 
sion (dashed line), and pure t'- J-dispersion t = (dotted 
line). This qualitative picture clearly demonstrates that 
due to the t'-terms the system becomes 'less strongly 
correlated' because of the 'liberating' of the bare car- 
rier. Figs. H and || show dispersions along the same lines. 
Fig. H describes the hole-doped system, and Fig. || the 
electron-doped one. Solid curves are related to the up- 
per and lower limit for the t' ,t" ,t N ,t s parameters from 
Table I, the dotted satellites are related to the same pa- 
rameters for the renormalized bare carrier weight (up to 
0.35). The dashed lines demonstrate the pure t-J model 
dispersions. 

Our calculations of the i'-part ( |l4] ) in the dispersion 
demonstrate that the 'bare' particle contribution for the 
Neel case is very close to the results for the Ising back- 
ground and that the 'dressed' particle contribution plays 
a minor role (< 0.2) compared to the 'bare' one. Using 
the data from Table I, the variations of the inverse mass 
values and the bandwidths for the hole and electron are 
found as follows: 



P'l = (3.8 - 5.3) J, 0\ = (3.5 - 5.0) J, 
/# = (1.3 -3.1) J, /3 2 e = -(0.7-1.6)J, 
W h = (5.1 - 7.1) J, W e — (5.2 - 8.4) J, 



(15) 



In spite of some changes in the inverse mass values (15), 
the features discussed earlier for the Ising case are still 
valid. Thus, the shifts of the band minima are rather 
large, the bandwidths are substantially wider as com- 



pared to those in the t-J model, and the role of the bare 
particle weight remains the most important. 

One of the unsolved problem of the t-J type mod- 
els is: whether the free-particle approximation work for 
the doped systems |3l| , |3^ , ^| . Therefore, the extension 
of the single-particle calculation onto the case of finite 
doping indeed requires justifications. Some of them can 
be found in Refs. |^,^l|,^8). We simply accept the free 
particle picture and have calculated the chemical poten- 
tial as the function of concentration for the dispersion 
law @, @. From Figs. |l|-| it is clear that the 
ground state is (±7r/2, ±7r/2) for the hole, and (0, ±7r), 
(±7r, 0) for the electron. Both minima are not shallow: 
|A(o,7r)-(7r/2,jr/2)| ~ J. At the same time, the bandwidths 
are larger than t-J ones (1 / F t_ ' 7 ~ 2J), and the accumu- 
lation of the DOS at low energy is smaller than for the 
t-J model case. This point may be of interest in view 
of intensive discussion of the possible Van Hove singular- 
ity at the optimal doping p5[ . The above mentioned 
two different tendencies to increase the minima shifts 
and to decrease the DOS require a quantitative consid- 
eration. Figs. [5] and || show our results for the DOS 
versus energy E/J and versus concentration 5, respec- 
tively. One can see that in spite of the DOS lowering 
it remains rather large at low energies (since (3 — 8) J 
<C 8t), and the peaks in the DOS shift to higher con- 
centrations (S m ax ~ 0.4 — 0.5), compared to the t-J case 
(S max ~ 0.1 - 0.2). 



IV. CONCLUSIONS 

To summarize, in this paper we have studied the mi- 
croscopically derived extended type of the t-t'-J model. 
Previous investigations of the three-band model allowed 
us to establish the necessary next-neighbor terms for the 
low-energy single-band model. In this work the efficient 
scheme, developed previously for the calculation of the t- 
J model parameters of the real Cu02 systems, has been 
applied to the effective t'-terms calculation. This has en- 
abled us to determine numerical ranges for all t'-terms 
both for the electron- and hole-doped systems. 

We also have performed simple calculations of the 
single-particle dispersion using the variational approach 
to the spin-polaron problem. Both types of doping have 
been considered within the above mentioned realistic val- 
ues of the t, t' e j f parameters. The importance of all t' e ^- 
terms has been clearly demonstrated: the t N , t terms of 
the second-order origin are responsible for the widening 
of the bandwidth, and the t', t"-terms arising from the 
0-0 hopping bring about the shift of the groundstate 
minima. It has been found that the minima shifts are 
large enough (~ J) and have opposite signs for the hole 
and electron systems. Also, the bandwidths have been 
found to be 1.5 — 4.0 times larger than those in the pure 
t-J model. 

An investigation of the finite doping regime for the 



5 



simple free-particle approximation has been carried out. 
The Van Hove peak in the DOS has been found to de- 
crease and move to a higher doping level from its t-J 
model position. 
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APPENDIX A: 

It will be useful for the future discussion to consider 
the simple limiting case of the three-band model, when 
the Cu-0 hopping is much less than both the Cu-Cu 
Coulomb repulsion \Jd and Cu-0 levels splitting A (t p d <C 
U d ,A). 

The additional hole moves over the spin background 
as a Zhang-Rice (ZR) singlet Q, whereas the electron 
moves as a vacancy. Corresponding hopping integrals in 
the lowest order of t p d/ A are: 



t 2 

t h = 0.53^ + 0.275i pp , t e 
t 2 

4 = 0.09-P - 0.124:t pp , t' e 
t 2 

t'l = 0.053-^ + 0.062i pp , t'l = - 



A 

A - 1.45t„ 



-(1.06 + 2.03^)4, 
A A 

_(0.l 8 -0.92^)V, 
V A ' A 

-(0.105 + 0.46^) V, 
A A 

(Al) 



The first two expressions for the hole hopping integrals 
coincide with those in the work by Jefferson et al [§, 
except the signs. The signs are opposite for the following 
reasons. 

In calculation of the transfer amplitudes of the t-J 
model singlet or vacancy from the three-band model, it 
is convenient to construct the basis of the wave func- 
tions as the direct product of the states at the sites: 



(A[|0)i <g) A||0) 2 ® • . .). The A\ operators can be con- 
sidered as the Hubbard operators at the site i, which 
create the singlet or spin state over the vacuum state 
|0)i. Namely, the two-site wave function of the form 
\s)i <8> | T)j has been used for the calculation of the singlet 
(\s)i) hopping integrals (see Refs. 0,||j28|). This conve- 
nience is due to the complex structure of the three-band 
low-energy states. For example, ZR singlet consists of 
the linear combination of a copper singlet (double copper 
occupancy), an oxygen one (double oxygen occupancy), 
and a copper-oxygen one. This state is hardly expressed 
in the terms of the Fermi creation operator which acts 
on the spin background. However, the use of the 'direct- 
product' basis is incorrect when more than one Hubbard 
A[ operators are of the Fermi type. This is evident from 
their anticommutativity. In the other words, the wave 
function of the half-filled background (one fermion per 
site) cannot be uniquely determined in this basis. For- 
mally, the two-site wave function |s)$(8>| T)j is determined 
correctly since the singlet |s)j is the boson. Neverthe- 
less, for the matrix elements of the hopping Hamiltonian, 
which consists of the creation and annihilation operators 
of an additional fermion, one gets: 



tij(1 



(kXT I), • (I T>^l) s . 



\s)i ® I T} 3 



MT |i®(T \i ■ I T)* ® I T) 3 



(A2) 



Thus, the matrix element reduces to the projection of the 
one half-filled state to the other and hence has an uncer- 
tain sign. Therefore, to avoid this uncertainty one has 
to use a more conventional basis, namely: c\\gs), where 
\gs) means the groundstate (one fermion per site), and c\ 
is the creation operator of the additional fermion. This 
careful approach leads to the signs of the hopping ampli- 
tudes as presented in Eq.( ]Al[ ) and in Table I. Note that 
this difference of the bases is absent for the single fermion 
in the lattice or for particles obeying Bose statistics. 

To obtain Eq. (Al) we have used the general expres- 
sions for tij from our previo us w orks JL2 28 1 . 

Approximate expressions (Al) for t, t' and t" are good 
enough for large and even moderate A. Thus, our calcu- 
lations of the t, t', t" realistic values for the Cu02 plane, 
show that 1 1'| < \t"\ at least for small values of V p d (Cu-0 
Coulomb repulsion). This is due to the partial compen- 
sation of the Cu-0 and 0-0 contribution to the i'-term, 
that can be seen in Eq.(Al). It should be noted that 



the t pp contribution to the hopping of the vacancy arises 
only in the third order of t p d/A (t pp /A). The above 
mentioned compensation of the 0-0 and Cu-0 ampli- 
tudes for the electron ig-term is even more pronounced 
than for the hole. The exact formulae fl2]| provide very 
small t' e for V p d = 0. 

To be more specific, in our calculations we have fol- 
lowed the idea of the narrowing of the uncertainty region 
for the low-energy model parameters using some experi- 
mentally observed values. This approach was developed 
in Ref. flq] and enabled us to calculate th and t e for the 
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real C11O2 planes. Our results for the t' ,t" ,t N ,t s hop- 
ping integrals obtained in the same way for the two types 
of doping are shown in Table I. 

We wish to stress the strong V^-dependence of *' both 
for the hole and electron carriers. This t'(V p d) depen- 
dence can be easily understood. The rising of V p d re- 
quires the lowering of the A for the fixed values of the 
observable quantities. Thus, the occupation number of 
oxygen sites also rises, which immediately leads to the 
growing of the t pp contribution. Note, that due to the 
smaller role of the oxygen degrees of freedom for the elec- 
tron, its hopping integrals are less varying than those for 
hole. This is in agreement with the cluster calculations 
by Eskes and Sawatzky 0,|2). In our three-band model 
calculations we used the experimental values of J = 0.14 
and J = 0.17 eV (lanthanum and yttrium systems). 



APPENDIX B: 



Explicit expressions for the components f k , Mk 01 are: 



v k = 



2S k 



-1 1/2 



Mk,<5 



xs k 

/47ke M 



(Bl) 



[YS k (A + 2S k )}^ 



■[(1 



(u + v)7 k e 



ikii 



where parameters X, Y, u, and v are expressed in terms 
of ground state correlators |l4| and for the Neel state 
are: A = 1.33J, X = 0.8, Y = 0.72, u = 0.42, v = 0.12. 
5 k = e*~ J (k) is the pure t-J model dispersion (13]). Since 
5 k ~ |*| at t > J, ^ k and // k , /ii k are weakly varying 
functions of k. 

Explicit expressions for the <5/3i ik , 8/32,k, ^k, -Bk coef- 
ficients (Eq.rtli)) are: 



A* = vlC\ + (/i k ) 2 Cf + {&) 2 Cl + ^ICl, 



(B2) 



where 



-5.8i" - hM N 
-8.2i"-8.2t 
9.8*' - Q.7t" 



N 



C\ = 3t' + 6t" + Ylt 

d 

Cf 

c{ 

C\ =-Zt' + 6t" 

C\ = 

d = 
d = 
c$ = 
cl = 



N 



■At* 



N 



12.8* 



■3t s 



-1.5t b , ci 
5.8*' + 
10.6*' - 
11.5*'- 
-5.8*' 



C 9 4 



N 



C} = 0, 



(B3) 



11.5*" + 23* N - 12. 2* 6 , 
- 16.4*" + 37.7* w - 15.1* 
- 23*" + 23*^ - 30.4* s , 
f 11.5*" + 4.6* s , 



-10.6*' + 16.4*" - 4.9*^ 
-11.5*' + 23*", 



3. Si- 



Numbers in Eq. @ result from (S* 4 Z S"|), {S?Sr) and 
other spin correlators for various neighbor sites. 
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Figure captions 



FIG. 4. Dispersion curves along the same lines as in Fig. 2 
for an electron. All notations as in Fig. 3. 

FIG. 5. DOS vs energy E/J for the pure t-J model and 
electron and hole t-t'-J model (upper and lower limits in Table 
I)- 



FIG. 6. DOS vs concentration 6. 



Tables 



Table I. t' , t" , t s ,t N hopping parameters for the hole and 
electron. First and second rows show the lower and upper 
limits, respectively. 





t'/\t\ 


t"/\t\ 


t N /\t\ 


t S /\A 


hole 


0.01 


0.12 


0.01 


-0.07 


t h /J = 2.5 


-0.25 


0.16 


0.07 


-0.16 


electron 


-0.09 


-0.11 


0.07 


-0.10 


t e /J = -3.0 


0.03 


-0.12 


0.11 


-0.12 



FIG. 1. Contour plot of e(k) for a hole. tu/J = 2.5, 
are taken as the average values of the upper and 
lower limit from Table I. Inner square is the magnetic Bril- 
louin zone boundary. r(0, 0), m[-k/2, tt/2), Z(ir, 7r), Y(0, 7r), 
and X(n, 0) points are indicated. 



FIG. 2. 
Dispersion curves along the lines T^M-^Z^Y^T, 
and Y — ► M — > X. Solid curve is the t-t'-J hole dispersion 
for the same parameters as in Fig. 1, dotted curve is the pure 
t'-J dispersion (t = 0), dashed curve is the pure t-J model 
dispersion. 



FIG. 3. Dispersion curves along the lines as in Fig. 2 for a 
hole. Solid curves correspond to the upper and lower limits of 
the ^-parameters (Table I). Dotted curves correspond to the 
same with renormalized bare carrier weight. Dashed curve is 
the pure t-J model dispersion. 
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